Geographically weighted generalized poisson regression model with the best kernel function in the case of the number of postpartum maternal mortality in east java

This study proposes a Geographically Weighted Generalized Poisson Regression (GWGPR) model with the best kernel function to obtain a model of the number of postpartum maternal mortality in East Java Province in 2020 and determine the factors that affect the number of maternal postpartum mortality in East Java in 2020. The kernel functions used in this study are fixed bisquare kernel, fixed tricube kernel, and adaptive bisquare kernel. Optimum bandwidth selection using the Cross-Validation (CV) method. The results obtained the best model is the GWGPR model with a fixed bisquare kernel because it produces the smallest AIC value of 194.92. Variables significantly affecting the number of maternal postpartum mortality in East Java in 2020 vary in each district/city where there are three regional groups. The percentage of pregnant women who had a pregnancy visit K1, the percentage of pregnant women who had a pregnancy visit K4, the percentage of households receiving cash assistance, and the ratio of hospitals and health centres have a significant effect on Kabupaten Blitar, Mojokerto, Gresik, Bangkalan, Blitar City, Mojokerto City, Surabaya City. While the five predictor variables together significantly affect districts/cities included in group 3, such as Ponorogo, Trenggalek, Tulungagung, Kediri, Malang, Lumajang, Jember, Banyuwangi, Bondowoso and so on. Some of the highlights of the proposed approach are:• Generalized Poisson regression model using the maximum likelihood estimation (MLE) method.• The kernel functions used in the Geographically Weighted Generalized Poisson Regression (GWGPR) model to determine bandwidth are fixed bisquare, fixed tricube, and adaptive bisquare kernel functions selected using the Cross Validation (CV) method.• The computation procedure is easy to implement.


Introduction
Regression analysis is a statistical method used to determine the relationship between variables. In a study, it is often found that the object of research is in the form of count data influenced by several explanatory variables. As the case of the number of maternal deaths is a rare event with discrete data, the appropriate regression analysis model to explain the relationship between the dependent variable in the form of discrete and rare data with independent variables in the form of discrete, continuous, categorical, or mixed data is the Poisson regression model. The Poisson regression model is a non-linear regression model used to model count or discrete data, with the main requirement being that the response variable has a Poisson distribution [1] . In the Poisson regression model, there is a specific assumption that the mean and variance of the response variable are equal (equidispersion).
However, there are times when the equidispersion assumption is only sometimes met. Namely, here can also be cases of overdispersion (the value of the data variance is greater than the value of the mean) or underdispersion (the value of the data variance is smaller than the value of the mean) in data modeling with a Poisson distribution [2] . One method that can be used to overcome these cases is Generalized Poisson Regression (GPR). The generalized regression model is one of the alternative models counting data with overdispersion or underdispersion cases in Poisson regression. Generalized Poisson Regression modeling produces a global regression model for all observation locations. In some instances, each observation location has a set of data that is different from one location to another, so to overcome this diversity, spatial data analysis can be used. Spatial regression modeling is performed on response variables in the form of count data that experience underdispersion or overdispersion and depend on the characteristics of the observed location using Geographically Weighted Generalized Poisson Regression (GWGPR). The GWGPR model develops generalized Poisson Regression by considering weights in the form of latitude and longitude of the observation location points. The GWGPR model's parameter estimator is local for each point or observation location. Some research on the GWGPR model is research by Ghanim, Al-Hasani, et al. (2021) entitled "Geographically weighted Poisson regression models with different kernels: application to road traffic accidents," analysing traffic accident data in the country of Oman using the Geographically Weighted Generalized Poisson Regression method with optimum weighting using adaptive kernel functions, namely box-car, bi-square, tricube, exponential and Gaussian [3] . Murakami, Daisuke et al. (2021) used the Geographically Weighted Generalized Poisson Regression method to analyze the covid-19 case in Tokyo, Japan. In this study, the bandwidth used was an adaptive and fixed Gaussian kernel [4] . Maternal Mortality Rate (MMR) is the number of women who die from a cause of death related to pregnancy disorders or their handling during pregnancy, childbirth, or in the maternal postpartum period (42 days after childbirth) but not caused by accidents or injuries [5] . MMR is an essential indicator in determining the degree of community welfare, especially in the health sector. MMR in developing countries such as Indonesia is still relatively high. Based on data from the Ministry of Health in 2020, maternal mortality in Indonesia reached 4627. This figure increased by 10.25% compared to 2019 [6] . East Java Province ranks second in the region with the highest MMR in Indonesia, with 565 maternal mortalities, most of which are caused by maternal postpartum mortality, amounting to 284 people [7] . In their research, Fronczak et al. (2005) stated that 75% of complications occurred postpartum, resulting in high maternal mortality [8] . This is an essential concern for the Indonesian government, especially in East Java Province, so MMR is one of the Sustainable Development Goals (SDGs) targets from 2015 to 2030, namely 70 per 100,000 live births in 2030 [9] . To reduce MMR, problems related to pregnancy, childbirth, and especially the postpartum period (42 days after delivery) cannot be separated from the various factors that influence it. This research will examine case data on maternal postpartum mortality in East Java in 2020 using the GWGPR method. In the GWGPR model, the optimum weight selection uses the fixed kernel bisquare, fixed kernel tricube and adaptive kernel bisquare functions, which are selected using the Cross-Validation (CV) method. Researchers use this method because there has been no previous research using the three kernels mentioned and the indicators considered to affect the dependent variable are also new.

Multicollinearity
Multicollinearity in the regression model can be determined using the Variance Inflation Factor (VIF), which is expressed in the following equation: 2 , SSR (Sum of Squares Regression) is the variation caused by the relationship between predictor variables, SST (Sum of Squares Total) is a measure of the variation of the value of X from the mean of X itself, X is the value of X at the th observation, X is average value of X and X is the predicted value of X at th observation.

Poisson distribution
In some research, it is often found that the object of research is in the form of count data influenced by several explanatory variables. A regression model based on the Poisson distribution can determine the relationship pattern between variables. The Poisson distribution is a discrete distribution that measures the probability of a certain number of events or events occurring within a certain period. The probability of a Poisson distribution with a random variable Y and parameter μ > 0 , which represents the number of successes that occur in a given time interval or region, is [10] where μ is the average number of 'successes' that occur in a given time interval or region and = 2 . 7183 . The Poisson distribution has the same mean value and variance as μ.

Poisson regression
The Poisson regression model is a non-linear regression model used to model count or discrete data, with the main requirement being that the response variable has a Poisson distribution. Poisson regression modeling the expected value of the response variable, with E(Y) as a linear function of the predictor variables ( ) is defined as follows [11] .
Since the mean (μ) of Poisson distribution is always positive, the function μ( X i , β) is chosen such that the linear predictor is: which can represent any real value to a positive real value. The logarithmic function is a suitable link function to model Poisson regression in the Generalized Linear Model approach, is: Thus, the Poisson regression model for the response variable Y, which is Poisson distributed with parameter μ is: Parameter estimation of the Poisson regression model can be done using the Maximum Likelihood Estimation (MLE) method. The maximum likelihood estimator for the parameter ( ) is expressed by ( ̂) , which is obtained by maximizing the likelihood function. Assuming that 1 , 2 ,…, are random variables that are mutually independent and ∼ Poisson ( μ ) , then the likelihood function for Poisson regression model is: The log-likelihood function of the Poisson regression model is: Furthermore, to obtain the estimator ( ̂) , the log-likelihood function in Eq. (9) ismaximized by lowering it concerning the parameter and equating it to zero. Parameter estimation using the MLE method obtained an equation that is not closed form, so we can estimate the parameter ̂b y using the Newton-Raphson iteration method as follows: where −1 ( ) is Hessian matrix and ( ) is gradient vector. Testing the parameters of the Poisson regression model is carried out to test whether the regression model parameters have a significant effect on the response variable (Y) . Model parameter testing using the partial and simultaneous tests is as follows. Hypothesis: The test statistic used: The L( ŵ ) function is the maximum likelihood value for a simple model with no predictors involved, and L ( Ω ) is the maximum likelihood value for a complete model with predictor variables involved. The rejection region of H 0 is to reject H 0 if the value of D( ̂) > 2 ( , ) , this means that there is at least one variable that has a significant influence on the response variable (Y) . The smaller the value of D( ̂) indicates, the smaller the error rate in the model [12] .
After testing the parameters simultaneously, we will continue with partial parameter testing. The partial parameter testing hypothesis is as follows: Hypothesis: The test statistic used: The rejection region of H 0 is to reject H 0 at the significance level if the value of |Z | > Z ( ∕2 ) this means that the variable p has a significant effect on the response variable.

Overdispersion or underdispersion in poisson regression models
The presence of overdispersion or underdispersion in Poisson regression can be detected using the devians, and Pearson Chi-Square values divided by the degrees of freedom. A value or quotient more significant than 1 indicates the presence of overdispersion ( var ( Y | ) > ( Y | ) ) , whereas if the result of dividing the two values is smaller than 1 indicating the presence of under dispersion var ( Y | ) > ( Y | ) , then one method that can be used to overcome this case is Generalized Poisson Regression [13] .

Generalized poisson regression
The Generalized Poisson Regression (GPR) model is used for counting data that experience equidispersion violations. The Generalized Poisson Regression model has the same form as Poisson Regression as follows: In the GPR model, the values of the parameters , β 0 , β 1 , β 2 , … , β will be estimated using the Maximum Likelihood Estimation (MLE) method with the likelihood The function of the GPR model are as follows: Next, Newton Raphson iteration is performed to maximize the log-likelihood function formulated as follows: Perform Newton Raphson iteration based on equation: The iteration process will stop if it has found an estimated value that converges to a value as follows: Furthermore, the estimation of the parameter using the Newton-Raphson iteration method is: A more accessible approach to estimating the parameter phi is to use the estimation of moments, equating the Pearson Chi-Square statistic with the degrees of freedom [14] . To determine with the method of moments can be expressed by the equation: Thus, the following iteration is obtained: where μ = exp ( ) > 0 GPR parameter testing is done with the Maximum Likelihood Ratio Test (MLRT) with the following hypothesis: The test statistic used: The H 0 rejection region is to reject H 0 if the value of D ( ̂) > 2 ( , ) , this means that there is at least one variable that has a significant effect on the model. If the simultaneous test decision is to reject H 0 , then the next step is to test the parameters partially to find out which parameters have a significant effect on the model. The hypothesis used is as follows: The test statistic used: The rejection region of H 0 is to reject H 0 at the significance level if the value of |W | > Z ( ∕2 ) this means that the variable p has a significant effect on the response variable.

Spatial heterogeneity test
Spatial heterogeneity occurs when the same predictor variable has an unequal effect on different locations within a study area. Spatial heterogeneity can be tested using the Breusch-Pagan (BP) test. Hypotheses used in the Breusch-Pagan test is as follows.

Geographically weighted regression
Geographically weighted regression is used to analyze spatial heterogeneity, where each parameter is calculated at each observation location, so each region or observation location has different regression parameter values. Systematically the Geographically Weighted Regression (GWR) model with response variable y and predictor variables 1 , 2 , … , at the th location can be written as follows: In spatial analysis, the role of bandwidth for GWR models is essential because the weight value represents the location of observations between one data and another. There are three types of bandwidth functions used in this study: The form of fixed bisquare is expressed by the following formula: ii. Fixed Tricube The following formula expresses the form of the fixed tricube: with ℎ is the bandwidth at the th observation location.
One of the methods used to determine the optimum bandwidth size is the Cross-Validation (CV) method. CV method is defined as follows: where Ŷ ≠ ( ℎ ) is the value of the estimator Y for observations at locations , not included in the calculation. The optimum bandwidth value ( ℎ ) is obtained when ℎ produces the minimum CV value.

Geographically weighted generalized poisson regression
Spatial regression modeling is performed on response variables in the form of count data that experience underdispersion or overdispersion and depend on the characteristics of the observed location using Geographically Weighted Generalized Poisson Regression. The difference between global Poisson regression and GWGPR is in parameter estimation. The global Poisson regression model has the same parameter estimator value for each observation location, so it is global. The GWGPR model produces a different parameter estimator value for each observation location area, so it is local [15] .
The GWGPR model is the development of Generalized Pois-son Regression by considering the weights in the form of latitude and longitude of the observation location points. The GWGPR model follows the distribution of Generalized Poisson Regression so that the probability function for each th location is as follows: GWGPR modeling uses a generalized linear model with a "g " function (link function) that connects the mean of the response variable with the linear predictor ( ). The link function used in the GWGPR model is the log link. The GWGPR model with as the latitude coordinate and as the longitude coordinate used as the parameter weight is

Parameter estimation of GWGPR model
In the GWGPR model, the method used to estimate the model parameters is the Maximum Likelihood Estimation (MLE) method, with the GWGPR model likelihood function as follows: Furthermore, the likelihood function in Eq. (28) is con-verted into natural logarithm form as follows: By substituting the value μ = exp X ( , ) the Eq. (30) is obtained.
)) ] The process of obtaining parameter estimators of the GWGPR model is by deriving the log-likelihood function for each parameter and then equating it to zero. Because the results are not close form, it is necessary to do a Newton-Raphson iteration with the following algorithm.

Determine the initial estimated values of the parameters
2. Form a gradient vector ( g ) with k estimated parameters 3. Form the Hessian matrix ( ) which is the second derivative of Eq. (30) 4. Substitute the values of ̂( , ) (0) into the elements of the vector and matrix to obtain the gradient vector ( ̂( , ) (0) ) and Hessian matrix ( ̂( , ) (0) ) . 5. Perform iteration using equation: 6. The iteration process starts from m = 0, and the value of ̂( u i , v i ) (m) is a set of parameter estimators that converge at the th iteration for the th location. 7. If you have not gotten a convergent parameter estimation, then the iteration continues back to step 5 until iteration = + 1 .
The iteration process is stopped if where is a very small number.

Simultaneous parameter significance test
The method used for simultaneous parameter testing in the GWGPR model is the Maximum Likelihood Ratio Test (MLRT) method. Testing the parameters of the GWGPR model is carried out to determine the significance of the β parameter with the following hypothesis.

Partial parameter significance test
Partial GWGPR model parameter testing is conducted to determine the effect of individual predictor variables on the response variable at each location. The partial parameter testing hypothesis is as follows.
The Statistics: where ( ̂ ( , ) ) is the standard error of β ( , ) . The H 0 rejection area is to reject H 0 at the α significance level if the value |Z | > Z ( ∕2 ) which means that variable has a significant effect on the response variable at each location in the GWGPR model.

Akaike information criterion
Akaike Information Criterion (AIC) is one of the criteria to determine the best model. A small AIC value indicates that the model is getting better. The best model selection using the AIC value criterion is as follows:

= −2 ln L ( ) + 2
is maximum log-likelihood model and is number of parameters estimated in the model.

Maternal postpartum mortality rate
Maternal Mortality Rate (MMR) is the number of female mortalities that occur during pregnancy or its handling during pregnancy, childbirth, or within 42 days after the end of preg nancy (postpartum period) but not mortality caused by accidents or injuries. MMR is an essential indicator in determining the welfare of society, especially in the health sector. The Ministry of Health in 2020 revealed that the number of maternal mortalities in Indonesia is still relatively high, reaching 4627 people. East Java Province ranks second with the number of maternal postpartum mortality at 284. The number of maternal postpartum mortality in East Java Province over five years is presented in the following figure.
Several studies have found that the factor that has a significant effect on the number of maternal postpartum mortality is the percentage of pregnant women who had a pregnancy visit K4 [16] . Furthermore, there is research on the number of maternal postpartum mortality, with one of the factors affecting maternal mortality is the percentage of handling obstetric complications in  each district/city [17] . Another study on the number of maternal postpartum mortality found that one of the factors affecting maternal postpartum mortality was the percentage of pregnant women who had a pregnancy visit K1 [18] . There is also research on factors that are thought to affect maternal mortality such as the percentage of households receiving cash assistance and the ratio of health centers and hospitals [19] . Based on the description above, in this study, several variables are used that are thought to affect the number of maternal postpartum mortality in East Java in 2020. Y is the number of maternal postpartum mortality, 1 is the percentage of pregnant women who had a pregnancy visit K1, 2 is the percentage of obstetric complications, 3 is the percentage of pregnant women who had a pregnancy visit K4, 4 is percentage of households receiving cash assistance, and 5 is ratio of health centers and hospital.

Characteristics of maternal postpartum mortality rate in east java
Descriptive analysis was carried out to obtain the characteristics of the number of maternal postpartum mortality (Y) and the factors that influence it. In 2020 the number of maternal postpartum mortality in East Java was 284. Descriptive statistics of all variables used and presented in Table 1 below. Table 1 . shows that the factors that are thought to affect the number of maternal postpartum mortality in East Java have a large enough variance and a large enough range that data heterogeneity is suspected. The predictor variable with the highest variance is the percentage of obstetric complications (  An overview of the distribution of cases of the number of maternal postpartum mortality in East Java in 2020 is presented in Fig. 2 .   Fig. 2. Maternal Postpartum Mortality Rate in East Java 2020.

Multicollinearity test
One of the requirements that must be met in Poisson regression analysis is multicollinearity detection. The presence of multicollinearity can be known from each predictor variable's VIF (Variance Inflation Factor) value. A VIF value of more than 10 indicates a case of multicollinearity. Table 2 shows that the VIF value for all predictor variables is less than ten, so there are no cases of multicollinearity in the data.

Distribution fit test
The distribution suitability test is used to determine whether the data is Poisson distributed or not. To conduct the test, you can use the Kolmogorov Smirnov Test with the following hypothesis: where F ( Y ) is the cumulative probability distribution. Test criteria: accept H 0 if the value of < D N , α in the Kolmogorov Smirnov table or the > α Based on the output, the − = 0 . 087 is obtained, so with a significance level of α = 0 . 05 , then the − = 0 . 087 > 0 . 05 , this means that the data on the number of maternal postpartum mortality is Poisson distributed.

Spatial heterogeneity test and the best bandwidth
Spatial heterogeneity testing is conducted to determine whether there are characteristics between observation location points. Spatial heterogeneity can be tested using the Breusch-Pagan (BP) test with the following hypothesis test: H 0 ∶ σ 2 1 = σ 2 2 = … = σ 2 38 = σ 2 (Variance between locations is the same) H 1 ∶ there is at least one σ 2 ≠ σ 2 ; = 1 , 2 , … , 38 (Variance between locations is the different) The test statistic value is 17.802 and the − is 0.003206. The significance level used is 5%, so 2 0 . 05;5 is obtained at 11.07. Therefore, it is concluded that reject H 0 , or the variance between locations is different. This means that there is spatial diversity between regions, or the characteristics between location points are different.
The existence of spatial effects causes spatial heterogeneity, so spatial weighting is needed. The best spatial weighting is obtained from the minimum Cross-Validation (CV) value's bandwidth value. Table 3 shows the optimum bandwidth selection value using the fixed bisquare, fixed tricube and adaptive bi-square kernel functions. Table 3 shows that the kernel function that produces the optimum bandwidth is the fixed bisquare Kernel function with a CV value of 1367.024. After obtaining the optimum bandwidth value, the spatial weight matrix for each observation location will be obtained by substituting the bandwidth value and Euclid distance.     Table 4 below. Table 4 shows that the AIC value of the Poisson regression model is 225. 19. The predictor variables that have a significant effect on the Poisson regression model are the percentage of handling obstetric complications ( 1 ) , the percentage of households receiving cash assistance ( 4 ) and the ratio of health centers and hospitals ( 5 ) . So that the effect of predictor variables 1 , 4 , 5 is significant but the effect of predictor variables 2 , 3 Table 4 is 72.026 with an independent degree of 32, resulting in a value of 2.25. This value is greater than 1, so the method used to overcome these cases is Generalized Poisson Regression (GPR). The Maximum Likelihood Estimation (MLE) method is used with Newton-Raphson iterations to obtain the estimated value of the GPR model parameters. Table 4 shows the estimated values of the GPR model parameters.
Based on Table 5 , it is known that the values of devians D ( ̂) is 202.4872. The significance level used is 0.05, so the value of 2 0 . 05;5 is 11.07. The devians value D ( ̂) is more than 2 0 . 05;5 which means reject H 0 . So it can be concluded that at least one predictor variable has a significant effect on the model.
The partial parameter test results show that the variables that have a significant effect are the percentage of households receiving cash assistance ( 4 ) and the ratio of health centers and hospitals ( 5 ) so that the GPR model for the case of the number of maternal mortalities in East Java Province in 2020 is as follows.

Geographically weighted generalized poisson regression
Testing the parameters of the GWGPR model is carried out to determine the significance of the parameters with the following hypothesis : Based on the result, the devians value D ( ̂) is 180.9181. The significance level used is 5%, so 2 0 . 05;5 is 11.07. The devians value of D ( ̂) is greater than the value of 2 0 . 05;5 so the test decision is to reject H 0 . This means that at least one predictor variable has a significant effect on the GWGPR model. Furthermore, the significance of the GWGPR model parameters is partially carried out to determine which predictor variables or factors have a significant effect in each region or location. The partial parameter testing hypothesis is as follows: The predictor variable is said to have a significant effect on the response variable if the value of |Z | > Z ( ∕2 ) . The significance level used is 5%, so the Z ( ∕2 ) value is 1.96. Based on the result, the percentage of pregnant women who had a pregnancy visit K1 ( 1 ) , the percentage of pregnant women who had a pregnancy visit K4 ( 3 ) and the percentage of households receiving cash assistance ( 4 ) had a significant effect on 38 districts/cities. While the percentage of obstetric complications ( 2 ) has a significant effect on 31 districts/cities and the ratio of health centers and hospitals ( 5 ) has a significant effect on 37 districts/cities. The significant variables in each district/city in East Java are presented as follows: • Variable 1 , 2 , 3 , 4 , 5 significant to 1 district (Kabupaten Pacitan). The grouping of districts/cities in East Java based on significant variables is presented in Fig. 3 . Based partial parameter testing, as an example, we will present parameter testing at the 3rd research location ( 3 , 3 ) = , namely Kabupaten Trenggalek, with parameter estimates shown in Table 6 . Table 6 shows that all predictor variables have a significant effect on the GWGPR model in Trenggalek District because the Z value is greater than Z 0 . 05∕2 . The GWGPR model formed in the case of the number of maternal postpartum mortality in the Kabupaten Trenggalek is as follows. Based on the GWGPR model in Kabupaten Trenggalek, every 5% increase in the percentage of pregnant women who had a pregnancy visit K1 will reduce the average number of maternal postpartum mortality in East Java by exp 0 . 00183 = 0 . 998 times, assuming other variables are constant. Suppose the percentage of pregnant women who had a pregnancy visit K4 increases by 5%. In that case, it will reduce the aver-age number of maternal postpartum mortality in East Java by exp ( −0 . 00583 ) = 0 . 994 times, assuming other variables are constant. The same interpretation applies to the predictor variables of the percentage of obstetric complications, the percent-age of households receiving cash assistance also variable the ratio of health centers and hospitals.
Determination of the best model. Based on Table 4 , it is known that for the Poisson regression model, three variables have a significant effect on the model, namely the percentage of obstetric complications ( 2 ) , the per-centage of households receiving cash assistance ( 4 ) also the ratio of hospitals and health centers ( 5 ) . Table 5 shows that the variables that significantly affect the GPR model are the percentage of households receiving cash assistance ( 4 ) and the ratio of hospitals and health, centers ( 5 ) . Meanwhile, in the GWGPR model, the variables of the percentage of pregnant women who had a pregnancy visit K1 ( 1 ) , the percentage of pregnant women who had a pregnancy visit k4 ( 3 ) , and the percentage of households receiving cash assistance ( 4 ) together influence 38 districts/cities. Further-more, the variables of obstetric complications and the ratio of hospitals and health centers have different influences for each district/city.
One of the criteria for determining the best model is to look at the Akaike Information Criterion (AIC) value obtained. The smaller the AIC value, the better the model used. The following are the AIC values for the Poisson, GPR and GWGPR regression models. Algorithm 1 Algorithm 1 is the procedure for estimating the parameters of the GWGPR model in the case of the number of maternal deaths in East Java in 2020. 1. Researchers are looking for case data on the number of maternal postpartum mortality in East Java Province in 2020 and longitude also latitude coordinate data. e. Calculate the best weight matrix using fixed bisquare kernel, tricube kernel and adaptive bisquare kernel.

Analyzing the
f. Parameter estimation of the GWGPR model using the Maximum Likelihood Estimation (MLE) method followed by Newton-Raphson iteration.
g. Simultaneous and partial significance testing of model parameters h. Interpreting the GWGPR model obtained.
i. Determining the best model using the AIC value, the best model is the model with the smallest AIC value.
5. The final result of this research is to get the best model to represent the number of maternal postpartum mortality and know the mapping of the number of maternal postpartum mortality in East Java Province.  The AIC value of the Poisson regression model is 225.19, the AIC value of the GPR model is 216.4872, and the GWGPR model with the best kernel (fixed bisquare) has an AIC value of 194.92. This means that the GWGPR model is the best because it has the smallest AIC value compared to other models ( Algorithm 1 , Fig. 1 , Table 7 ).

Conclusion
The average number of maternal postpartum mortality in each district/city in East Java Province was 7.47 in 2020, with the highest death of 30 cases in Kabupaten Jember and the lowest case in Sampang, Kediri city, Madiun City and Batu city. The high variance value of 31.553 indicates a reasonably high difference in maternal postpartum mortality in each district/city. The best kernel function is the fixed bisqure kernel function with a CV value of 1367.024. Based on the results obtained, the GWGPR model is the best model to model the case of the number of maternal postpartum mortality in East Java in 2020 because it produces the smallest AIC value of 194.92 with three groupings obtained based on the results of partial testing for GWGPR modeling.

Ethics statements
The data used in this study are secondary data derived from the 2020 East Java Provincial Health Profile published by the East Java Provincial Health Office and some secondary data derived from the official website of the Badan Pusat Statistik (BPS) East Java.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data Availability
Data will be made available on request.